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THE TYPE IIB SUPERNOVA 201 1DH FROM A SUPERGIANT PROGENITOR 
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ABSTRACT 

A set of hydrodynamical models based on stellar evolutionary progenitors is used to study the nature 
of SN 201 ldh. Our modeling suggests that a large progenitor star — with R ~ 200 Rq — , is needed to 
reproduce the early light curve of SN 201 ldh. This is consistent with the suggestion that the yellow 
super-giant star detected at the location of the SN in deep pre-explosion images is the progenitor star. 
From the main peak of the bolomctric light curve and expansion velocities we constrain the mass of 
the ejecta to be ~ 2 M , the explosion energy to be E = 6 — 10 x 10 50 erg, and the 56 Ni mass to be 
approximately 0.06 Mq. The progenitor star was composed of a helium core of 3 to 4 Mq and a thin 
hydrogen-rich envelope of ~ 0.1 Mq with a main sequence mass estimated to be in the range of 12-15 
Mq. Our models rule out progenitors with helium-core masses larger than 8 Mq, which correspond 
to Mzams ^ 25 Mq. This suggests that a single star evolutionary scenario for SN 201 ldh is unlikely. 
Subject headings: hydrodynamics — supernovae: general — supernovae: individual: SN 201 ldh 



1. INTRODUCTION 

Type lib supernovae (SNe lib) are transitional ob- 
jects within the family of core-collapse SNe (CCSNe), 
as their spectroscopic classification evolves from Type II 
(i.e. with H lines), to type lb (i.e. dominated by helium 
lines). SNe lib were rec ognized as a new class with the 
disc overies of SN 1987K (lFilippenkoiri 988') and SN 1993J 
(e.g. iNomoto et al.l 119931 : lFilippenkolll997t ). As the rest 
of CCSNe — which includes H-rich types II-P and II-L, H- 
less type lb, and He-less type Ic — , they are believed to 
arise from the violent death of stars with initial masses 
greater than 8 Mq. Massive stars may suffer consid- 
erable mass loss during their evolution, due to strong 
stellar winds or mass transfer to a binary companion. 
Therefore, the vast spectroscopic and photometric diver- 
sity observed among CCSNe is related to the ability of 
the progenitor star to retain its outermost layers before 
the explosion. Despite the efforts to improve our under- 
standing of the progenitor of each subclass of CCSNe, 
many questions remain open. In this sense, the study of 
a very well-observed object can shed light on the nature 
of CCSNe and their massive progenitor systems. 

SN 201 ldh was discovered on May 31.893 by am- 
ateur astronomers and immediately confirmed by the 
Palomar Transi ent Factory (PTF) in the nearby spi- 
ral galaxy M51 (iGriga et al.ll20111 [Silverman et alJi20lH 
lArcavi et all 120111 ) which had also hosted three other 
CCSNe in the past 17 years. Strong constraints on 
the date of explosion, to better than 0.6 days were es- 
tablis hed using pre- and post-SN imaging (|Arcavi et al.l 
12011 . SN 2011dh was extensively monitored in a wide 
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wavelen gth range, including ea rly detections in radio and 
X-r avs (ISoderberg et al.ll2012D . It was classified as type 
lib {A rcavi et al.ll2011h . a relatively rare subclass of CC- 
SNe, based on the optical spectrum. Soon after dis- 
covery, a source was identified as a possible progeni- 
tor of SN 2011dh in archival, multi-ban d HST images 
(jMaund et al.ll2Qlll IVan Dvk et ai] 120111 ). Photometry 
of the source is compatible with a yellow super-giant 
(YSG) star. 

The detection of pre-supernova (pre-SN) objects in 
high-resolution imaging has provided important infor- 
mation on the possible progenitors of several SNe. For 
SNe II-P, these observations confirm the red supergiant 
nature of the progenitor, as previously suggested by 
theory. An upper limit of Mzam s = 16.5 ± 1.5 Mq 
was sugg ested for this type of SNe (jSmartt et alJ [2009). 
Ho wever, such upper lim i t has been recently revised 
bv iWalmswell fc Eldridgi (j2012D . They found a higher 
value for the maximum mass of SN II-P progenitors of 
Mzams = 21 Mq when additional extinction due to 
dust produced in the red supergiant wind is taken into 
account. For other subtypes of CCSNe, detections have 
not been as common and only in a few cases have the pro- 
genitors be en conclusively ident ified, e.g the type II-pec 
SN 1987A dGilmozzi et alll987lh th e type lib SN 1993J 
(|Aldering et al.|ll994t iMaund et all 120041) . and the type 
fin SN 2005gl ({Gal- Yam et alJl2007f ).~ 

To determine the main-sequence mass (Mzams) from 
pre-SN imaging, it is necessary to derive a luminosity and 
intrinsic colour from the photometry, and compare with 
so me evolutionary trac k. T his was done for SN 2 011dh 
bv lMaund et~aT1 poTTI ) and lVan Dvk et al.l ([2011b . But 
although both studies obtained consistent effective tern- 
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perature and luminosity, and they employed the same 
evolutionary model, the derived Mzams was different. 
Because color uncertainties are expected to arise from un- 
known mass-loss history of the progenitor, iMaund et al.1 
(| 2 1 If) assumed only the luminosity was reliable, and de- 
rived Mzams = 13 ±3 M from the end point of the evo- 
lutionary track that m atched the luminosity. Meanwhile 
iVan Dvk et~aT1 (j2011ft derived M ZA ms = 17-19 M by 
choosing the closest track that matches the luminosity 
and color of the source in the Hertzsprung- Russell (H-R) 
diagram although this point does not correspond with the 
final positi on of the star at the end of its evolution. Al- 
ternatively, iMurphv et al.l (|2011f) compared stellar popu- 
lation synthesis with the stellar association surrounding 
SN 2011dh to derive Mzams- Assuming that the stars 
in the vicinity of the SN are coeval, they concluded that 
th e value of Mzams is most likely close to the estimate 
of IMaund et all (12ml . 

Some authors have suggested that the YSG star de- 
tected in the pre-SN images is not the actual progenitor 
of SN 201 ldh but its binary companion or even an un- 
related object. These authors suggest, instead, that the 
exploding object was a compact star. The arguments for 
this are based, first, on the shock veloc ity derived from 
radio and sub-millimeter observations (ISoderberg et al.l 
I20H iBietenholz et~all [20121: iKrauss et al.l 12012ft . The 
relatively high shock velocity found, v sw ~ 0.1c, would 
indicate a compact pro genitor, or a type cllb SN, i n 
the scheme proposed by IChevalier fc Soderbersj (|2010f ). 
The second argument for a comp act progenitor was in- 
troduced by lArcavi et al.l (|2011ft and is based on the 
quick decline of the optical light curve (LC) soon af- 
ter the explosion, as compared with SN 1993J, along 
with a low temperature derived from an early-time 
spectrum, as compared w ith analytic expressions by 
IRabinak fc Waxmanl f|201lD . Such a compact progenitor 
would be inconsiste nt with the YSG star ide ntified by 
IMaund et all (|2011l ) and lVan Dvk et al.l (|2011l ). which is 
expected to have a radius of « 270 i? (from log L — 4.92 
Lq, T eff = 6000 if, given L = i? 2 T c 4 ff in solar units). 

Early observations such as those available for 
SN 201 ldh provide a unique opportunity to analyze the 
physical properties of the progenitor. In particular, one 
of the most direct ways to estimate the size of the pro- 
genitor is by modeling the LC during the cooling phase 
that occurs after the shock break-out and before the re- 
heating by radioactive decay. Observations during this 
phase are very scarce due to its short duration and so 
they are very valuable. The studies described above do 
not present a specific modeling of the early light curve. 
With the aim of assessing the nature of the progenitor 
we set out to perform a more detailed modeling of the 
available observations. 

A well-known method to estimate the properties of the 
progenitor object, as well as the explosion energy, and 
the amount and distribution of the radioactive material, 
is to compare the observations with predictions from hy- 
drodynamical models. In this paper, we present a set 
of hydrodynamical models applied to stellar evolution- 
ary progenitors that aim to elucidate the compact or ex- 
tended nature of the progenitor of SN 201 ldh, as well 
as the main physical parameters of the explosion. Our 
initial and LC models are presented in § [2] A compari- 



son between models and observations is done in § The 
global physical parameters are studied in § 13.11 and the 
sensitivity of the LCs on the initial radius is investigated 
in § 13.21 In §[4] we discuss the results and summarize our 
main conclusions § [5] 

2. STELLAR AND SUPERNOVA MODELS 

To study the nature of SN 201 ldh we compare the 
observations with a set of hydrodynamical models ap- 
plied to initial structures from stellar evolutionary calcu- 
lations. Our supernova models are computed using a one- 
dimensional Lagrangian h ydrodynamic code w ith flux- 
limited radiation diffusion (Ber sten et al.|[20Tl[) . The ex- 
plosion itself is simulated by injecting a certain amount 
of energy near the center of the progenitor object, which 
produces a powerful shock wave that propagates through 
the progenitor transforming thermal and kinetic energy 
of the matter into energy that can be radiated from the 
stellar surface. During the propagation of the shock 
wave, explosive nucleosynthesis produces unstable iso- 
topes of iron-group elements, mainly 56 Ni. The decay 
of 56 Ni — > 56 Co — >- 56 Fe produces energetic 7-rays and 
positrons that are thermalized, providing extra energy 
that contributes to power the LC. The code includes 7- 
ray transfer in gray approximation for any distribution 
of 56 Ni , assuming a constant value for the 7-ray opacity, 
« 7 = 0.06 y e cm 2 g~ 1 , where y f . is th e number of elec- 
trons per baryon ([Swartz et al.l 11995ft . Note that this 
approximation is appropri ate for typica l SN c onditions 
as it has been shown by iSwartz et al.l (|1995| ) through 
comparison with Monte Carlo simulations. Nucleosyn- 
thesis is, however, not consistently calculated but it is 
included as a pre-explosion condition, assuming that the 
energy released does not contribute significantly to the 
dynamics of the explosion itse lf, and that it only effect is 
on the chemi cal structure (see lWooslev fc Weaver lfl990l ; 
lArnettl 119961 for further details of this assumption). 

LCs of SNe lib have been successfully reproduced us- 
ing helium stars with a very t hin hydrogen envelope 
(< lMd) as pre-SN models dShigevama et al.1 119941 : 
iWooslev et all 119941 : iBlinnikov et all 119981 ) . Therefore, 
we adopt the same type of model as input for our hydro- 
dynamic cah2u]ations^_Tlie_Fle core m odels were calcu- 
lated bv lNomoto fc Hashimoto I (|1988l ) using a single-star 
evolutionary code that follows the evolution of the He 
core from He burning until the collapse of the core assum- 
ing solar initial abundance. The external H-rich enve- 
lope of the He core was replaced with the proper bound- 
ary conditions obtained from the hydrostatic and ther- 
mal equilibrium enve lope models (jNomoto &; Sugimoto I 
119721: iNomoto IH974T ). To take into account the thin hy- 
drogen envelope, we smoothly attached a low-mass H- 
rich envelope in hydrostatic and thermal equilibrium to 
the He core. The presence of the external envelope sig- 
nificantly modifies the radius but not the mass of the 
progenitor. 

Specifically, we employ four different initial models 
with He core masses of 3.3 M Q (He3.3), 4 M (Hc4), 
5 Mq (He5), and 8 M Q (He8), which correspond to 
the stellar evolution of single stars with main-sequence 
masses of 12, 15, 18, and 25 Mp, respectiv ely (see 
iSugimoto fc Nomoto 1 119801: iTanaka et all 120091 for the 
He core mass - main-sequence mass relation) . Figure Q] 
shows the initial density profile as a function of radius 
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for the He4 model and for models with different H-rich 
envelopes attached to the core of the He4 model. The 
latter models will be analyzed in § 13.21 

The chemical structure of these He stars before the ex- 
plosion consists of a Fe core of ~ 1.4 M Q , surrounded 
by concentric layers of Si-rich, O-rich, and He-rich mate- 
rial. For each of the initial models, explosive nucleosyn- 
thesis was calculated with a r eaction network that in- 
clude s 280 isotopes up to 79 Br (|Hix fe Thielemannlll996l 
1999). The results of the nucleosynthesis depend on the 
progenitor mass and the kinetic energy of the explosion 
and they have already be en presented in previous studies 
(| Thielemann et all 1 19961) . An important characteristic 
of the models is the variation in the mass of the O-rich 
layer which goes from 0.2 M Q for the He3.3 model to 
3M Q for the He8 model. Therefore, the determination 
of the oxygen mass using late-time spectra may provide 
a constraint on the mai n-sequence mass of the progenitor 
(|Kawabata et al.ll2010l) . 

In Section 13.11 we study the explosion models for the 
three least massive initial models (He3.3, He4, and Hc5), 
using several values of explosion energy, 56 Ni mass, and 
56 Ni distribution. In this analysis the mixing of 56 Ni 
is modified from that of the nucleosynthesis calcula- 
tions in order to explain the rising part of the observed 
LC, as fre quently done to mode l the LC of SNe Ibc 
(e.g. see iShigevama. et"ai1 119941: iWooslev et all 119941 : 
iBlinnikov et all 120001 ) . Such a mixing between the core 
and the He layer in Type lb and lib supernovae is caused 
by th e Ray leigh- Taylor ins t ability at the base of the H e 
layer ijHachisu et al. II1991L fl99i llwamoto et al. Ill997l ). 
Actually, more recent 3D numerical calculations predict 
even larger mixing of 56 Ni than previous 2 D calculations 
(|Hammer et al.ll2010Hjoggerst et al.ll2010D . This result is 
in better agreement with the extensive mixing observed 
in e.g. SN 1987A. In Section [3J2] we study the effect of 
the external H-rich envelope which must be present to 
explain the H spectral features observed at early times. 
Finally, the He8 model is discussed in Section |4] to test 
the plausibility of a massive single progenitor star. 

3. LIGHT CURVE MODELING 

The g'-band LC of SN 2011dh showed a quick decline 
during 1-3 days after an initial peak with M^ cak — —16.5 
mag. This initial decline was followed by a re-brightcning 
in all bands leading to a second peak with bolometric lu- 
minosity of L bo i ~ 1-4 x 10 42 erg s -1 (Af bo i = -16.6 
mag) at t « 20 days. Double-peaked LC have been 
observed for very few SNe, one of which was the well- 
known SN 1993J, although with a slower initial decline 
and larger luminosity, as compared with SN 201 ldh. The 
scarcity of observations during the early declining phase 
is due to its very short duration, which makes it diffi- 
cult to catch. Theoretically, however, all CCSNe should 
present this early behavior which is a consequence of 
cooling after the shock breaks out from the surface of 
the object. The duration of the cooling phase is mainly 
regulated by the size of the progenitor, although there 
is also a dependence on energy, mass and envelope com- 
position. More compact structures, such as Wolf-Rayet 
stars, produce faster declines than extended progenitors, 
such as red supergiants. 

On the other hand, the heat source that powers the 



second peak is provided by the radioactive decay of 56 Ni 
produced during the explosive nucleosynthesis and its 
daughter 56 Co. Modeling the LC around the second 
peak provides information about the explosion energy 
(-E), ejecta mass (M e j), and the mass and distribution of 
56 Ni . Such global properties for SN 2011dh are analyzed 
in Section 13.11 while the effect of the progenitor radius 
which essentially affects the early decline phase of the 
LC is studied in Section [3T2T 

The observed bolometric LC of SN 2011dh and the 
expansion velocities for the Fe I I lines used in this pa- 
per are t aken from lErgon et alj (|2012l ). A distance of 
7.1 Mpc (|Takats fe Vinkol 120061) . a Galactic reddening 
of E(B — V) = 0-035 mag and no host-galaxy redden- 
ing (|Arcavi "et~aD[20Tlh were assumed in the calculations. 
For the early cooling phase, there are only a few observa- 
tions available in a limited wavelength range, which pre- 
vents us from calculating a bolometric luminosity during 
that epoch. Hence we dir e ctly a dopt the g'-band data 
published by lArcavi et al.l ()2011[ ) . In the present work 
we assume the explosion time was 2011 May 31.5 UT, 
that is the mid-point between the last non-detection and 
the discovery date. 

3.1. Global properties of SN 201 ldh 

To estimate physical parameters, such as explosion en- 
ergy (E), ejected mass (M C j), and the mass ( 56 Nimass) 
and distribution of 56 Ni, we analyze the bolometric LC 
around the second peak, and the photospheric velocitiy 
i v ph) evolution. For this purpose, the He star models 
presented in section [2] are adopted as initial configura- 
tions. To save computation time, no external H envelope 
was added at this point because this does not affect the 
LC around the main peak. Therefore, the progenitor is 
assumed to be compact with a radius of R ss Rq. 

First we analyze the general dependence of the LC and 
u p h on variations of the physical pa rameters for mode l 
He4. This is the He core used bv lShigevama et "all ()1994[ ) 
to model the bolometric LC of SN 1993 J with E m 1.2 
foe (1 foe= 1 x 10 51 erg), 56 Ni mass of ~0.08 M , and 
assuming substantial mixing of 56 Ni out to the He lay- 
ers (see also our own model of SN 1993J in Figure [T51 
Section[0]). Following that work, for SN 2011dh we test 
physical parameters near those values. 

Figures [2] and [3] show the effect of the explosion energy 
on the LC and on v p h, respectively, for energies of E = 
0.8, 1.0, and 1.5 foe. More energetic explosions produce 
larger kinetic and radiative energy, which is reflected in 
higher photospheric velocity and global brightness. After 
the cooling phase, all LCs evolve through a minimum 
with luminosity L m i n , followed by a broad second peak 
with maximum luminosity L pe ak, and then they enter 
a more or less linear decline, or "tail" phase. In the 
following analysis we denote the epochs of L m - ln and £ pea k 
as < m in and £ pca k, respectively. From Figure [5] we see 
that L m - m and L poa k increase with E while the luminosity 
during the tail phase shows a dependence in the opposite 
direction. Note that t m - m is more or less independent of 
explosion energy, but £ pea k decreases with E. Finally, 
models with lower E produce broader second peaks. 

The effect of the 56 Ni mass on the LC is shown in Fig- 
ure@for 56 Ni mass = 0.05, 0.06, and 0.07 M . The LCs 
are remarkably equal until about 10 days before i pca k 
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when the models begin to differ in the sense that larger 
56 Ni mass produces larger luminosity and wider LC. 
£peak and the tail luminosity increase with 56 Ni mass 
but no appreciable effect is seen in L min . Also note that 
t mm and i pea k are essentially insensitive to the amount 
of synthesized 56 Ni. 

The distribution of 56 Ni could be different than the 
one predicted by one-dimensional nucleosynthesis calcu- 
lations where it is concentrated in the inner regions of the 
progenitor. Mixing of 56 Ni out to the helium envelope 
is expected to occur because of Rayleigh- Taylor instabil- 
ities during the propagation of the shock wave which are 
not properly taken into account in our one-dimensional 
prescription. To study the effect of mixing, we calculated 
three models with 56 Ni linearly mixed (in mass coordi- 
nate) out to 75% (He4Mix75), 85% (He4Mix85), and 95% 
(He4Mix95) of the total initial mass. Note that the dis- 
tributions of 56 Ni adopted are in good co ncordance with 
recen t 3D calculations (see e.g. Fig. 6 of Hammer et al.l 
l2010h . 

The resulting LCs are shown in Figure [5] The main 
effect of mixing is seen at t < i pca k although there is also 
a smaller effect on the tail luminosity. For more extensive 
mixing, t min is produced at earlier times because of the 
earlier heating by radioactive decay. This also produces 
a higher L m ; n , while the value of -L pea k seems not to be 
very sensitive to mixing. Instead the effect is more visible 
for the LC width which increases with more extended 
mixing. 

The results of these tests show that the effect on v p h 
of the 56 Ni mass and distribution is essentially negligi- 
ble. Therefore, the kinetic energy available in the ejecta 
is the only relevant parameter that determines the pho- 
tospheric velocity evolution, which imposes strict con- 
straints on E for a given M c - y 

This comparative analysis allows us to find a set of pa- 
rameters that reproduce the observations of SN 201 ldh. 
For He4, these parameters are E = 1 foe, and 
56 Nimass= 0.065 M , mixed out to 95% of the total 
mass or out to 9000 km s _1 in velocity space. Note that 
although 56 Ni was mixed far out in the ejecta, since the 
adopted mixing function was linear, the amount of 56 Ni 
in the outer layers was small (< 1.3 x 10 -3 Mq for v > 
5000 km s _1 ). We adopted a cut mass (M cut ) of 1.5 Mq 
that is assumed to form a compact remnant. Therefore, 
the mass of the ejecta is M e j = M to tai — -Mcut = 2.5M Q . 
A summary of the parameters used in this and the follow- 
ing calculations is presented in Tabic [TJ where we show 
the main-sequence mass, He core mass, and radius of the 
progenitor, the mass cut assumed, the ejecta mass, the 
explosion energy and the 56 Ni mass. 

We also tested the He3.3 and He5 models with smaller 
and larger pre-SN mass, respectively. The parameters 
found for these cases are slightly different from those of 
the He4 model (see Tabled]), although the 56 Ni mixing 
was kept up to 95% of the total mass. For instance, in 
order to keep the agreement with the observed velocities, 
we modified the explosion energy so that the ratio E/M C j 
remains approximately constant. Figures |5] and [7] show 
the LCs and v p h, respectively, for the parameters of mod- 
els He3.3, He4, and He5 compared with the observations. 
From the figures, it is clear that the three models give 
good fits to the observations, with the two least massive 



ones producing slightly better results. 

Models He3.3 and He4 give the best fits to velocity and 
LC, respectively. This means that there is a slight tension 
in the solution of the progenitor properties based on the 
comparison of our calculations with these observables. 
However, given the uncertainties involved, it is likely that 
the best global fit arises from a model with parameters 
in the range of these two models. That is, M e j = 1.8-2.5 
M Q , E = 0.6-1.0 foe, and 56 Ni mass = 0.060 -0.065 Mq 
mixed out to a velocity of w 9000 km s _1 . Such a progen- 
itor with a He core mass between 3 and 4 Mq corresponds 
to a star of ~ 12-15 Mq on the main sequence. Although 
there are uncertanties in the determination of the main 
sequence mass from the He core mass, depending for ex- 
ample on convective overshooting, rotation, etc., for the 
range of He core masses considered here different stel- 
lar calculations provide similar results (see Figure 1 of 
Tanaka et al. 2009). Stars of such initial masses have 
no known mechanism to almost completely remove the 
H envelope, as required for producing a SN lib, except 
through binary interaction. In Section [4] we study the 
possibility of more massive progenitors and discuss the 
binary scenario. 

3.2. Progenitor radius 

To test the effect of the progenitor radius on the LC, 
it is critical to obtain observations during the adiabatic 
cooling phase before the radioactive decay becomes the 
main source to power the LC. For the case of SN 2011dh, 
only a few data points in the g'-band are available at 
those times. Here we compare these data with two hydro- 
dynamical models: a compact one with an initial radius 
of ks 2 Rq, and an extended one with R « 270 Rq, consis- 
tent with the L and T e g of the YSG star detected in pre- 
SN images. Both models have a He core mass of 4 Mq 
and the same physical parameters as model He4 (see Ta- 
ble [J) . The extended model (He4R270) was constructed 
by smoothly attaching to the helium core a H/He en- 
velope in hydrostatic and thermal equilibrium (see Fig- 
ure [1]). The mass (M cnv ) and the helium mass fraction 
(Y) of the envelope for a uniform composition were fit so 
that the L and T e g became equal to the estimate from 
pre-SN imaging. This yielded values of M cnv = 0.1 Mq 
and Y — 0.8. Note that the total H mass in this model 
is ps 0.02 M<7), similar to the value derived from spectral 
modeling ( Arcavi et al.ll201lh . 

Figurc[8] (left panel) shows the bolometric LCs for both 
the compact and the extended progenitor models. The 
effect of radius is only noticeable before (k5 days, which 
is roughly where the bolometric LC data begin. At t ~ 1 
day the difference in bolometric luminosity between the 
two models is very large, of up to one order of magnitude. 
This difference is mostly due to the extra amount of en- 
ergy required to expand a more compact structure. For- 
tunately, earlier data obtained in the g' band are avail- 
able to discriminate between the two models. Therefore, 
to compare with the g'-band data we calculated synthetic 
LCs assuming black-body emission and integrating the 
flux through the corresponding filter transmission func- 
tion. Figure [8] (right panel) shows a comparison of our 
models with the g'-band observations for compact (blue 
line) and extended (red line) progenitors. Again, both 
models differ significantly in the g' band at t < 5 days. 
While the extended model produces a pronounced spike, 
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in agreement with the observations, the compact progen- 
itor shows a much weaker bump. From this analysis we 
conclude that an extended progenitor is favored by the 
early-time observations. 

Further evidence for this scenario can be sought by 
analyzing the temperature of the ejecta. Figure [9] shows 
the evolution of effective temperature for both models. 
There are no important differences in the effective tem- 
per i ature_for_£> 2 days , contrary to what was suggested 
bv lArcavi et al.l fj201 lh . In that paper, the authors es- 
timated a black-body temperature of « 7600 K from 
the spectrum of SN 2011dh at 2.8 days, and compared 
it with the effective temperature derived from an an- 
alytic expression of (|Rabinak fc Waxmanl 120 111 RW11 
hereafter) that strongly depends on the progenitor ra- 
dius, and less strongly on other parameters. Assuming 
a radius of ^10 13 cm at t = 2.8 days after the explo- 
sion, the predicted temperature is T e s » 17000 K. This 
discrepancy is one of the main arguments they used to 
rule out an extended progenitor, such as a YSG. On the 
contrary, our numerical calculations give values of T e g = 
7730 K and 7690 K at this epoch for compact and ex- 
tended progenitor, respectively, which is compatible with 
the value derived from the spectrum. 

An additional model named He4R270NH is shown in 
Figure IH1 (dotted line). This model is similar to He4R270, 
but with larger helium mass fraction in the envelope 
of Y w 1, i.e. essentially hydrogen- free. One can 
see that the presence of H in the envelope of the ex- 
tended progenitor (model He4R270) improves the agree- 
ment with the observation by reducing the tempera- 
ture, as compared with a model where no H is included 
(model He4R270NH). The composition of the envelope 
thus plays an important role in the determination of the 
effective temperature. 

For comparison, Figure |9] also shows the behavior of 
T e g calculated using equation (13) of RW11 for the cases 
of compact and extended progenitors. The adopted radii 
are the same as those employed for the hydrodynamical 
calculations, i.e. R = 2AR Q and R = 270 Rq. Note that 
only at very early times do the analytic models agree 
with the hydrodynamical calculations. At later times 
the agreement breaks down. This could be partially due 
to the effect of recombination in the ejecta which is not 
considered in expression (13) of RW11. Another possible 
reason for the disagreement could be differences in the 
initial density structures between both formulations (see 
Figure Q] for our initial density profile). A more detailed 
study is needed to clarify the origin of the discrepancy, 
but that goes beyond the scope of the present work. In 
any case, it is important to bear in mind that the analytic 
expressions were applied at an epoch (« 2.8 days) that 
could be outside their validity range. 

In any case, the T e e is not directly comparable with 
the black-body temperature derived from the spec- 
trum. A more direct comparison can be done using the 
color temperature (Tg). Following the prescription of 
lEnsman fc Burrows! (|1992[ ) we estimated Tq as the tem- 
perature at the "thermalization" deptlQ, which led to 
values of 8500 K and 8300 K at 2.8 days for models He4 

The "thermalization" depth is calculated as the layer where 
3 r a bs T sct ~ 1 , where r sc t is the optical depth for scattering and 
■Tabs is the optical depth for absorption. 



and He4R270, respectively. Although these values are 
somewhat higher than the value estimated from the spec- 
trum, the discrepancy is not important given the uncer- 
tainties in the time of explosion (~ 0.6 day) and in the 
estimations of the color temperatures. Because of the 
small differences in temperature found at t ~ 2 days be- 
tween compact and expended progenitors, the available 
temperature measurement is not a suitable discriminator 
between these scenarios. 

Finally, we analyzed whether it is possible to improve 
the comparison between models and early observations 
assuming different values of the progenitor radius than 
that inferred for the YSG star. Figure [TU] shows the 
bolometric (left panel) and (/'-band (right panel) LCs for 
models with progenitor radii of 50, 100, 150 and 200 
Rq. All of these configurations have the same He core 
taken from the He4 model, and they were constructed in 
a similar way as He4R270, i.e. smoothly attaching an Fi- 
nch envelope to the core (see Figured]). We denote these 
models as He4R50, He4R100, He4R150 and He4R200. As 
seen from the figure, it is clear that models with R ps 200 
Rq are more consistent with the early-time data. This 
finding is not affected by the systematic uncertainty in 
the luminosity that would arise from an error of 1 Mpc 
in the distance. 

We conclude this analysis by claiming that a progenitor 
with radius similar to that of a YSG star, as suggested 
from pre-SN detections, is compatible with the early ob- 
servations of SN 2011dh. Moreover, we find that radii 
much smaller than 200 Rq fail to reproduce the obser- 
vations. 

4. DISCUSSION 
4.1. Single versus binary progenitor 

SNe lib require the hydrogen-rich envelope of the pro- 
genitor star to be almost completely removed before the 
explosion. Two alternative mechanisms of envelope re- 
moval have been proposed to explain the progenitors of 
SNe lib, lb, and Ic, thereby called "stripped-envelope 
SNe": (1) strong stellar winds in massive single stars, 
and (2) mass transfer in close binary systems. In the first 
scenario a very massive star with a main sequence mass 
> 30 M q is required for the mass-loss rate to be large 
enough (jHeeer et al.l [200l iGeorgy et al.l 120091) . This 
type of star has a He core mass > 8 Mq previous to the 
explosion. The upper limit of the main-sequence mass 
may be even l arger according to recent stellar wind mass- 
loss rates (see IBouret et al]l2005t lEldridge fc Vinkl 12005 
iFullerton et al.ll2006f) . Iiithe binary scenario, less mas- 
sive stars are allowed with He core masses prior to the 
explo s ion in the range o f 3-6 Mq (jPodsiadlowski et"afl 
[19931: lYoon et all l201Ch . In the previous section we 
showed that such He-core mass range is in very good 
agreement with the observations of SN 201 ldh. 

To further test the possibility of a single-star progen- 
itor, we calculated a model based on a progenitor with 
main sequence mass of 25 Mq which forms a He core of 
8 Mq previous to the explosion (we call this model Hc8). 
In Figures [6] and [7] we show the LC and v p h, respectively, 
for model He8 using the same 56 Ni mass and distribu- 
tion as found for the He4 model of section 13.11 but with 
a larger explosion energy of E = 2 foe in order to repro- 
duce the peak luminosity. Clearly, this model does not 
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agree well with the observations. While decreasing the 
explosion energy can improve the match to the expansion 
velocities, it would worsen the fit to the LC irrespective 
of the 56 Ni mixing assumed. Note that the timing of 
the second peak imposes an important constraint on the 
He core mass. More massive helium stars reach the LC 
maximum at later times because the heat produced by 
radioactive decays takes longer to diffuse out. The He8 
model is too massive to produce the second maximum at 
« 20 days as observed for SN 2011dh, even assuming the 
most extreme 56 Ni mixing. 

This situation cannot be remedied by assuming a dif- 
ferent distance to M51 within the uncertainty of 1 Mpc. 
This would move all the LC points systematically up- 
ward or downward by 0.12 dex, but it would not change 
the shape of the LC, i.e. its width and peak timing. 
Therefore our hydrodynamical modeling indicates that 
the He core mass must be lower than 8 M Q , which favors 
a binary origin for SN 2011dh. 

Recently it has been suggested that the binary chan- 
nel is preferred to explain most of the stripped-envelope 
SNe. This is supported for example by po pulation studies 
comp ared with the observed rates of SNe ( Eldridae et al.l 
[20081 : ISmith et all l20"TIh . In the case of SNe lib the 
binary picture is further strengthened for two reasons: 
Firstly, stellar evolution of single stars requires a precise 
fine tuning of the initial parameters to leave a thin H en- 
velope previous to the explosion. T his is much more nat- 
urally explained in a binary context fl P o dsiadlowsk i et al.1 
[19931 IWooslev et al.l 11991 lYoon et alJl201Ct ). Secondly, 
the detection of the binary companion in the case of 
the type lib SN 1993J, and possibly also in the case 
of SN 2001ig. A K0 super giant star was identified as 
the progenitor of SN 1993 J (|Aldering et al.lll994f >. How- 
ever, the pre-SN photometry showed an excess in the UV 
and B bands that was associated with a blue companion 
star. Approximately a decade after the explosion the 
blue supergiant companion w as confirmed (|Maund et al.l 
[200llMaund fc Smartdl20"o"9h . In the case o f SN 2001ig 
a poss ible companion star was reported by iRvder et al.1 
(2006]). 

The observed pre-SN photometry of SN 201 ldh fits 
well the spec tral energy distribution (SEP) of a sin - 
gle YSG star (IMaund et al.ll20ll IVanDvk et al.ll20Tl . 
Therefore, for the binary sceneario to apply, the com- 
panion flux should have no appreciable effect on the pre- 
SN photometry. We have performed binary evolution 
calculations to test the possibility of systems that are 
compatible with the pre-SN observations. For this pur- 
po se, we employed the int e ractin g-binary code developed 
bv lBenvenuto fc De Vitol ((200 3?). The code was initially 
constructed for low-mass systems and recently extended 
to compute the evolution of massive stars by including 
nuclear reactions until the end of oxygen burning. Fur- 
ther details of t hese calculations are pr esented in our 
companion paper iBenvenuto et al I (|2012ft . 

Figure [TT] shows the evolutionary tracks in the H-R 
diagram for both components of a binary system com- 
posed of a primary (donor) star with an initial mass 
of 16 Mq, a secondary (accreting) star of 10 M , and 
an initial orbital period of 150 days. The conservative 
mass-transfer case (0/3 = 1) is shown in the left panel of 

2 /3 represents the fraction of material lost by the donor star that 



the figure, where all the matter lost by the primary star 
is accreted onto the secondary. The right panel shows 
a non-conservative case with an accretion efficiency of 
25% (/3 — 0.25). With this configuration, independently 
of the adopted value of /3, the primary star ends its evo- 
lution within the region of the H-R diagram defined by 
the pre-SN photometry and its uncertainty. Also note 
that the mass of the primary star prior to the explosion 
is « 4M Q , which is consistent with the hydrodynamical 
modeling presented in the previous section. 

Furthermore, in all the above calculations some hy- 
drogen remains in the envelope of the primary star at 
the moment of the explosion with a similar mass of 
» 4 x 1O~ 3 M . This H mass is high enough to be de- 
tected in the spectra an d therefore it is cons istent with 
the SN lib classification (jDessart et al.ll2011| ). 

To study the effect of the secondary on the pre- 
explosion phot ometry we used atmosphere models of 
iKurucz I |l993)0 for the corresponding T e g, surface grav- 
ity and radius of both binary components to compute 
synthetic magnitudes through the observed HST band- 
passes. As can be seen in Figure [TU the secondary star 
is significantly bluer than the primary at the time of ex- 
plosion. Thus, the largest contribution to the flux of the 
system from the secondary is in the bluest band, F336W. 
For both values of /3 the contribution of the secondary 
in the rest of the bands is < 4%. In the conservative 
case (P — 1) the secondary would increase the flux of 
the system in the F336W band by 50%, producing a 
marginal 1.5 a detection, considering the measurement 
uncertainty in this band. In the non-conservative case 
with j3 = 0.25 the contribution of the secondary in the 
bluest band is only of 30% of the total flux which falls 
well inside the photometric uncertainty (0.6 a). 

If the binary scenario really applies to the case of 
SN 201 ldh, the companion star should be recovered in 
future observations, when the SN becomes faint enough. 
From the synthetic photometry describ ed above and as- 
sumi ng a distance to M51 of 7.1 Mpc (jTakats &; Vinkol 
120061 ) and no host-galaxy extinction, the expected bright- 
ness of the secondary in the F336W band is 24.9 mag and 
24.2 mag for /3 = 0.25 and (3 = 1, respectively. In the 
redder bands, the brightness lies within 26-27 mag. A 
search for the secondary star using the HST when the 
SN fades enough is feasible and can serve to ultimately 
test the binary scenario. 

4.2. Compact versus extended progenitor 

In the literature, a more compact progenitor has been 
suggested for SN 20 11dh than the YSG object detected 
at the SN position (lArcavi et al.l [201 It iSoderberg et al.l 
[2011 I Van Dvk et al.l 120 111 ). This was based on (1) a 
comparison between the early decline rate in the LC of 
SN 1993 J and SN 201 ldh, (2) an incompatibility between 
the temperature derived from an early-time spectrum 
and that estimated from analytic expressions for an ex- 
tended object, and (3) large shock velocity derived from 
radio observations. 

Regarding the first point, while it is true that SN 1993J 
showed a slower decline during the cooling phase as corn- 
is accreted by the secondary. 

3 The Kurucz stellar atmospheres atlas was downloaded from 
http: / /www. stsci . edu/hst/observatory/cdbs/k93models .html 
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pared with SN 201 ldh, this can be explained by the dif- 
ference in size of the proposed progenitors. SN 1993J was 
suggested to have a red supergiant pro genitor with a ra- 
dius of « 600i? Q (|Maund et al.ll200l logL/L© = 5.1 
and logT c ff = 3.63), whereas the YSG star that was pro- 
posed as the progenitor of SN 201 ldh has a smaller ra- 
dius of ps 200i?©. To test this statement we calculated a 
model for SN 1993J based on the He4 progenitor with an 
envelope attached to make L and T c g consi stent with the 
pre-e xplosion photometry of SN 1993 J fsee lMaund et all 
[2001 . Fi gure [12] shows a compari son between our mode l 
and the observations of SN 1993J (jRichmond et al.lll994t l 
assuming the explosion date was 1993 March 27.5 UT. 
An explosion energy of E = 1.2 foe and a 56 Ni mass 
of 0.084 M© were used in this simulat ion. Such values 
are co nsistent with those adopted by IShigcvam a et al.1 
<|1994D . We also show the He4R270 model and the ob- 
servations of SN 201 ldh for a clear comparison between 
both SNe and models. From the figure it is clear that the 
difference in radius is enough to explain the difference in 
decline rate of the early light curves between both SNe. 

With respect to the temperature discrepancy of the 
second point above, using hydrodynamical models in 
§ 13.21 we did not find any incompatibility with the spec- 
trum temperature at t » 2.8 days when assuming an 
extended object. Note that our models show that the dif- 
ferences in temperature for t > 2 days between compact 
and extended progenitors are too small to discriminate 
between these scenarios. In our modeling we are able to 
determine the radius and thus set a lower limit of about 
200 Rq (see Section 13. 2[) based on the comparison with 
the early LC and not with the temperature. 

Regarding the issue of the radio properties, the 
Shockwave velocity appears to be too high for an 
extended object based on the picture presented by 
IChevalier fc Soderbergl (|2010| ). However, the connection 
between the shock velocity and the physical radius of 
the proge nitor is not direct and r el ies on several assump- 
tions (see ISoderberg et~aT1 ((20121 ): IKrauss et al.1 (I2012T) ). 
Therefore, we think that this argument alone cannot be 
used to rule out the extended progenitor scenario. More- 
over, it is worth noting that the mass-loss rate e stimated 
from radio observations by IKrauss et al.l (|2012f) is com- 
patible with a YSG progenitor and that no variability in 
the radio light curve, as has been observed for some com - 
pact SNe lib (jRvder et alJ 120041 : ISoderberg et al]l200fih . 
was reported to date. 

5. CONCLUSIONS 

We have calculated a set of hydrodynamical models 
applied to stellar-evolution progenitors in order to study 
the nature of SN 2011dh. Comparing our models with 
the observed bolometric LC during the second peak and 
with line expansion velocities we found that a progenitor 
with He core mass of 3.3-4 Mq, an explosion energy of 
6 - 10 x 10 50 erg, and a 56 Ni mass of « 0.06 M Q repro- 
duce very well the observations assuming a distance of 
7.1 Mpc to M51. This type of model is consistent with 
a main sequence mass between 12 and 15 Mq. Remark- 
ably, the range of mass found with our hydrodynamical 
modeling is in very good agreement with the estimates 
from two other independent methods, i.e. pre-SN imag- 
ing and stellar population analysis. This is different from 



the situation generally encountered for SNe IIP where 
LC modeling estimates main sequence masses that are 
higher up to a f actor of two than those estimated from 
pre-SN imaging (jUtrobin fc Chugaill200l lSmarttll2009l ). 

We have studied the effect of the progenitor radius 
on the early LC and temperature evolution. We found 
that a progenitor with radius similar to that of the YSG 
star detected in the pre-SN images is compatible with 
the early observations of SN 201 ldh without contradic- 
tion with the temperature that is derived from the spec- 
trum, as opposed to what lArcavi et al.1 fj201 lh found us- 
ing analytic models. Furthermore, progenitors with radii 
< 200 Rq fail to reproduce the early g'-band LC. Al- 
though our hydrodynamical models show differences in 
the temperature evolution between extended and com- 
pact progenitors, these differences are less marked than 
those predicted by analytic expressions and they are al- 
most unnoticeable for t > 2 days. Therefore, the spec- 
trum temperature at t « 2.8 days is not useful in this 
case to discriminate between compact and extended pro- 
genitors. 

We have tested and ruled out progenitors with He core 
masses > 8M© which correspond to Mzams ^ 25 Mq. 
Considering the limitations at such stellar masses of 
single-star winds to expel the H-rich envelope almost en- 
tirely, as required for SNe lib, this result is in favor of a 
binary origin for SN 201 ldh. 

We have also performed binary evolution calculations 
with mass transfer to test the possibility of systems 
that are compatible with the pre-SN observations of 
SN 2011dh. We have shown that a system with 16 Mq + 
10 Mq and an initial period of 150 days predicts that the 
primary star ends its evolution in the H-R diagram at the 
right position as com pared with the YSG star detected 
in the pre-SN images (jMaund et al.l[20Tll [Van Dv k et ahl 
120111 ). Furthermore, the He mass of the primary at the 
end of the evolution was ~ 4Af Q , which is consistent 
with our hydrodynamical modeling. The binary evo- 
lution calculations further predict that some hydrogen 
mass of ~4x 1O _3 M is left in the envelope of the pri- 
mary star, which is required to produce a SN lib. 

To test the binary scenario we have studied the ef- 
fect of the putative companion star on the pre-explosion 
photometry in comparison with the observations. Be- 
cause the secondary star is predicted to be much hotter 
than the primary star, we found that the largest effect 
appears in the blue and UV filters. The contribution of 
the secondary to the flux in the F336W band, however, is 
marginal, at the 1.5 a level. The contribution is further 
decreased to the 0.6 a level when non-conservative mass 
accretion is considered. However, our predictions can be 
tested in a few years time by a search for a blue object 
at the location of the SN. 
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Fig. 1. — The initial density distributions as a function of radius for model He4 (red line) and for models with different H-rich envelopes 
attached to the core of He4. 



10 



Bersten et al. 




Fig. 2. — Sensitivity of the bolometric LC to changes in the explosion energy. The He4 initial model (see § 13. Il l and three different values 
of the explosion energy, E = 0.8, 1.0, 1.5 foe, were used in these simulations. The observed bolometric LC of SN 2011dh (points) is shown 
for comparison. 
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Fig. 3. — Sensitivity of the photospheric velocity evolution on the explosion energy. The He4 initial model (see § 13,11) and three different 
values of the explosion energy, E = 0.8, 1, 1.5 foe, were used in these simulations. Fc fl line velocities measured from spectra of SN 2011dh 
(points) are shown for comparison. 



12 



Bersten et al. 



42.5 




20 40 60 80 

t [days] 

Fig. 4. — Sensitivity of the bolometric LC to changes in 56 Ni mass. The He4 initial model (see § 13. Il l and three different values of the 
56 Ni mass, 56 Nimass = 0.05,0.06,0.07 Mq, were used in these simulations. The observed bolometric LC of SN 2011dh (points) is shown 
for comparison. 



SN 2011dh 



13 



42.5 



~i 1 r 



n 1 r 



~i 1 r 



~i 1 r 



42 



w 

W 

CD 



o 41.5 



He4Mix95 
He4Mix85 
He4Mix75 




41 



_i i i_ 







_l I L_ 



20 



_i i i_ 



40 
t [days] 



_1 I L_ 



60 



80 



Fig , 5. — Sensitivity of the bolometric LC to changes in the 56 Ni distribution. The He4 initial model with 56 Nimass= 0.06 Mq (see 
§ETT} and three different degrees of mixing, up to 75% (He4Mix75), 85% (He4Mix85), and 95% (He4Mix95) of the total initial mass, were 
used in these simulations. The observed bolometric LC of SN 201 ldh (points) is shown for comparison. 
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Fig. 6. — Observed bolometric LC of SN 2011dh (points) compared with the results of the LC calculations for models He3.3 (magenta 
solid line), He4 (red solid line) and He5 (blue solid line). An extra model, He8 (green solid line) is also included to show that larger helium 
core mass is not compatible with the observations. The error bar indicates the size of the systematic uncertainty in luminosity as derived 
from an uncertainty of 1 Mpc in the distance to M51. The physical parameters used in each simulation are given in Table[l] 



SN 2011dh 



15 




Fig. 7. — Evolution of the photospheric velocity for models He3.3 (magenta solid line), He4 (red solid line), He5 (blue solid line), and 
He8 (green solid line) compared with measured Fe II line velocities of SN 2011dh. The physical parameters used in each simulation are 
given in Table 1. Note that model with larger helium mass overestimate the observed velocities. 
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Fig. 8. — Observ ed and modeled bo lometric LCs (left panel) and absolute (/-ba nd LCs (right panel). The dots show the observed 
bolometric LC from Ergon ct al. (2012). and the (/-band LC from Arcavi ct al. (2011). The blue solid lines show the results for the compact 
progenitor model He4. The red solid lines corresponds to an extended progenitor, He4R270, consistent with a YSG star as detected in 
pre-SN photometry of SN 201 ldh. Clearly, the effect of the progenitor radius is only important before the radioactive material becomes 
the main source of radiative power. The larger radius is necessary to reproduce the early part of the (/-band LC. The error bars indicate 
the size of the systematic uncertainty that corresponds to an uncertainty of 1 Mpc in the distance to M51. 
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Fig. 9. — Time evolution of the effective temperature for the compact model He4 and the extended mod els He4R270 and He4R270NH 
(see § 13.21 1, The effective temperature calculated using the analytic expression of Rabinak & Waxman (2011) with R = 270 Rq (solid black 
line) and R = 2.4 Rq (dashed black line) are also shown. The black-body temperature (cyan dot) estimated from a spectrum of SN 201 ldh 
obtained at 2.8 days is included for comparison. 
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Fig. 10. — Bolometric LCs (left panel) and (/-band LCs (rig ht panel) for mo dels with the same explosion energy a s our preferred 
model, but different initial radii. The observed bolometric LC IIErgon et alj |20121 and g'-band LC l lArcavi et alj|201lf ) of SN 2011dh 
(cyan dots) are shown for comparison in each panel. The error bars indicate the size of the systematic uncertainty that corresponds to an 
uncertainty of 1 Mpc in the distance to M51. The radius variation is accomplished by attaching essentially massless (< 0.01 Mg) envelopes 
to the He4 model. Larger radii produce higher early luminosity for t < 5 days but no appreciable effect is seen at later times. 




Fig. 11. — Evolutionary tracks in the Hertzsprung- Russell diagram for both components of a binary system with solar composition. The 
stars have Mzams of 16 Mq and 10 Mq and an initial period of 150 days. Left panel: assuming conservative mass transfer (/3 = 1). 
Right panel: non-conservative mass transfer (/3 = 0.25). Labeled dots along the tracks indicate the masses of the stars (in solar masses) 
while mass-transfer by Roche-Lobe overflow occurs. The primary star (black line) ends its evolution with a mass ~ 4 Mq and with effective 
temperature and luminosity consistent with the YSG star detected in pre-SN images (red rectangle), independently of the value of /3. 
However, the evolution of the secondary star (blue line) strongly depends on the assumed accretion efficiency. In the conservative case, 
the bolometric luminosity of both stars at the moment of the explosion of the primary star (big dot) is similar, while for /3 = 0.25 the 
secondary is ft: 2 mag weaker than the primary at the end point. 
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Fig. 12. — Co mparison between theoretical bolometric LCs (solid lines) and observations (dots) for SN 1993J (Richmond ct al. 1994) 
and SN 2011dh (Ergon ct al. 2012). In both cases a progenitor was used with He core mass of 4 Mq. Envelopes of different radii were 
attached to the core to make L and T e ff consistent with pre-explosion imaging of each SN. The differences during the cooling phase are 
well explained as differences in the size of the progenitor while the differences in the main peak and tail are related with different different 
values of the explosion energy and 56 Ni mass. 
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TABLE 1 

Physical parameters of the explosion models 



Model 


M ms a 
[Mq] 


[Mq] 


R* 
[«©] 


M cut 
[Mq] 


Mcj 
[Mq] 


E 
[foe] 


56 Nimass 
[Mq] 


Hc3.3 


12 


3.3 


2.5 


1.5 


1.8 


0.6 


0.065 


He4 


15 


4. 


2.4 


1.5 


2.5 


1. 


0.06 


He5 


18 


5. 


1.8 


1.6 


3.4 


1.2 


0.065 


He8 


25 


8. 


1.3 


1.8 


6.2 


2. 


0.065 



a Main sequence mass 



